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We study the conductivity of a class of disordered continuum systems represented by the Swiss- 
cheese model, where the conducting medium is the space between randomly placed spherical holes, 
near the percolation threshold. This model can be mapped onto a bond percolation model where 
the conductance a of randomly occupied bonds is drawn from a probability distribution of the 
form a"". Employing the methods of renormalized field theory we show to arbitrary order in e- 
expansion that the critical conductivity exponent of the Swiss-cheese model is given by t^'~^{a) = 
(d — 2)u + max[(^, (1 — a)~^], where d is the spatial dimension and u and (ji denote the critical 
exponents for the percolation correlation length and resistance, respectively. Our result confirms a 
conjecture which is based on the 'nodes, links, and blobs' picture of percolation clusters. 

PACS numbers; 64.60.Ak, 05.60.-k, 72.80.Ng 
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I. INTRODUCTION 

Percolation [Q is one of the best studied problems in 
statistical physics, both because of its fundamental na- 
ture and its vast array of applications. The most natural 
type of percolation, perhaps, is continuum percolation, 
where the positions of the constituting elements are not 
restricted to the discrete sites or bonds of a regular lat- 
tice. A simple example of continuum percolation is a 
conducting material with uniformly-sized holes placed at 
random. Due to its similarity to Swiss cheese, this model 
is commonly called the Swiss-cheese model. 

Since the holes are allowed to overlap, the system 
ceases to support electrical transport when the total vol- 
ume of the holes exceeds a critical fraction Qc- Near this 
percolation threshold qc the conducting network consists 
of many narrow bottlenecks each of which is bounded 
by inter-penetrating holes. Thus, it is plausible that the 
Swiss-cheese model can be mapped onto the random re- 
sistor network (RRN) problem where conducting nearest 
neighbor bonds on a hyper-cubic d-dimensional lattice 
are randomly occupied with a probability p. Apparently, 
the bottlenecks are playing a role similar to the occupied 
conducting bonds. However, the bottlenecks have a wide 
distribution of widths, in contrast to the standard RRN, 
where all occupied bonds are identical. Due to the wide 
distribution of neck-widths the Swiss-cheese model corre- 
sponds to a modified RRN in which the conductances a of 
the individual occupied bonds have a broad distribution 
in the form of a power law cr"" with < a < 1 1^. Due 
to its relation to the Swiss-cheese model, we abbreviate 
such a RRN by RRN^^. 

It is well established that the purely geometrical per- 
colation exponents for the Swiss-cheese model are in con- 
formity with their analogs in the discrete models . For 
example, the correlation length ^ is governed in both 
models by the same exponent v. The reason is that only 



the connectivity of the bottlenecks is relevant for the geo- 
metrical exponents. As the widths of the individual bot- 
tlenecks do not matter, they can be regarded identical in 
the context of connectivity properties and the problem is 
essentially equivalent to standard discrete percolation. 

The situation is different for critical exponents pertain- 
ing to transport quantities. Let us consider the conduc- 
tivity exponent t for the RRN. It describes the decrease of 
the average macroscopic conductivity E when the critical 
occupation probability pc is approached from above M] 



(1.1) 



The conductivity exponent is related to the resistance 
exponent (j) governing the average resistance Mn between 



two terminal sites 
cluster M, p| 



and x' known to be on the same 



Mr{x,x') ~ \x-x'\ 
via the scaling relation 

{d-2)iy + 4>. 



(1.2) 



(1.3) 



The conductivity exponent for the RRN^*~^ on the other 
hand depends on a, i.e.. 



^sc 



(1.4) 



Early estimates of t^'-'(a) were given by Kogut and Stra- 
ley [7], and Ben-Mizrahi and Bergman [||. Later Stra- 
argued based on the 'nodes, links, and blobs' pic- 
10[ of percolation clusters that t^'^{a) is given by 



ley I 
turc 



f^{a) = (d - 2)v + max [0, (1 - a)' 



(1.5) 



(see also Machta et al. |T^). Without relying on the as- 
sumptions of the 'nodes, links, and blobs' picture t^'^{a) 
has been addressed by Lubensky and Tremblay (LT) llj] 
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from a renormalization group (RG) perspective. After 
some controversy their perturbation calculation to 
first order in the deviation from the upper critical di- 
mension for percolation e = 6 — d shows agreement with 
Eq. (|L5l). 

The paper in hand presents our field theoretic study 
of the conductivity of the RRN^^ . Our analysis builds 
up on the field theoretic RRN^^ Hamiltonian by LT. We 
discuss the RG flow for the whole regime < a < 1 to 
arbitrary order in a diagrammatic expansion. The central 
result of our work is that Eq. (1.5) holds to arbitrary 
order in e. 

The outline of the remainder of this paper is the fol- 
lowing. Section || describes the modeling. We define the 
percolation problem under consideration. Then we show 
how the average resistance and the related average con- 
ductance can be derived from a generating function. We 
explain how the replica trick facilitates averaging and 
leads to an effective Hamiltonian. Next, this effective 
Hamiltonian is refined into a field theoretic functional. A 
scaling analysis concludes Sec. ||. It reveals the relevance 
of the field theoretic couplings associated with the con- 
ductances of the occupied bonds. In Sections III and IV 



we actually compute the generating function for the aver- 
age conductance by employing field theory augmented by 
renormalization. A Gell-Mann-Low RG equation (RGE) 
provides us with the scaling behavior of the average con- 
ductance near criticality. Our analysis is partitioned into 
two cases. In Sec. Ill we consider a — and basically 



review the known results for the RRN. Section IV deals 
with the case of prime interest, viz. < a < 1. Finally, 
concluding remarks are given in Section [V. Technical 
details are relegated to Appendices ^ andjB. 



Since we are going to calculate the conductivity expo- 
nent i^^(a) via the average conductance M^^i we need a 
precise definition of this quantity. Commonly this defini- 
tion is based on a setup in which a fixed external current 
/ is applied between two leads at lattice sites x and x' 
known to be on the same cluster. In this setup one could 
measure the resistance R{x,x') between the two termi- 
nals and then average with respect to g, 



Ml['ix,x')^{R{x,x% 
with the average being defined as 



(••■>: 



{xix,x')---)g 

{x{x,x'))g 



(2.2) 



(2.3) 



Here x{x,x') is an indicator function that takes on the 
value one if x and x' are located on the same cluster 
and zero otherwise. The average macroscopic resistance 
( p.2[ ), however, has the severe drawback that it is not well 
defined, because the average bond resistance diverges. 
Thus, it is preferable to work with the average conduc- 
tance instead which is given by 



M|C,(x,x') = (i?(x,a;')~')' 



(2.4) 



Note that x probes only geometrical connectivity. 
Hence, {x{x,x'))g can be identified with {x{x,x'))c, 
where {■■■)c denotes averaging over all diluted lattice 
configurations C of the corresponding standard bond per- 
colation model. Accordingly {xix, x'))g is nothing more 
than the usual percolation correlation function, i.e., the 
probability P{x,x') that x and x' are connected. 



II. THE MODEL 

We are about to consider a bond percolation model 
on a d-dimensional hyper-cubic lattice where the con- 
ductances of the occupied bonds are independently and 
identically distributed random variables. To be specific, 
the distribution function g of the conductance Cf, of any 
bond b is taken to be 

0(a) = (l-p)5(a)+pf(a) (2.1a) 

where 

f(a)=.(l-a)ao-i(^^) " (2.1b) 

with (7 G [0, (To] and < a < 1. For the relation of 
the RRN^^ defined by this choice to continuum percola- 
tion we refer to Halperin et al. [Q. Note that f has the 
important feature that the average resistance of an occu- 
pied bond is infinite |lj] . This is a key distinction to the 
standard RRN and also to a RRN with noise modeled 
by a narrow distribution of bond conductances (cf., e.g., 
Ref. and references therein). 



A. Generating function 

In this section we review how one can devise a generat- 
ing function for M^5i based on the ideas of Stephen [ p^ . 
We demonstrate that this generating function indeed 
serves its purpose and explain how the average conduc- 
tance can be extracted from it. 

Stephen introduced the quantity 

^3j(x) =exp(iA-T4) , 0. (2.5) 

Vx — {Vx'^\ ■ ■ ■ ,Vx^^) is a _D-fold replicated variant of 
the voltage Vx at lattice site x and A ~ (A'^^\ ■ • • , A'^-') 
is, apart from a factor — i, a replicated external current. 
The corresponding scalar product is defined as A ■ T4 = 
'}2^=i V^i^^A^"^ The physical content of ip^ix) will be 
explained below. 

In order to proceed towards the desired generating 
function we now consider the two-point correlation func- 
tion of tp^{x) 

g(x,x',a) =(v3j(x)^_x(x'))^^ (2.6) 
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where the average is defined by 



orthogonaHty relations 



7-D 



D 



(a) 



X exp 



j a = l 



(2.7) 



The product over j is taken over all lattice sites. Z is a 
normalization factor given by 



j exp 



-lPi{V}) 



P{{V}) denotes the dissipated power 

p{{v}) = E '^b^b = E '^-^ - 

b <i.j> 



(2.8) 



(2.9) 



with the summations running over all bonds. Vb abbrevi- 
ates Vi—Vj, where i and j are the lattice sites belonging to 
the respective bond, and accordingly, aij = Ub- P{{V}) 
is the replicated version of the power with all voltages 
replaced by their replica ted analogs. 

Before evaluating Eq. (2.6) we need to comment on reg- 
ularization issues. First, it is imp orta nt to realize that 
the integrands in Eqs. (2.7) and (2.8) depend only on 
voltage differences and hence the integrals are divergent. 
To give these integrals a well defined meaning one can 
introduce an additional power term Physi- 
cally the new term corresponds to grounding each lattice 
site by a capacitor of unit capacity. The original situ- 
ation may be restored by taking the limit of vanishing 
frequency, w ^ 0. Second, it is not guaranteed that 
Z stays finite because infinite voltage drops may occur. 
Thus, the limit lim/j^o is not well defined. This prob- 
lem can be regularized by switching to voltage variables 9 
taking discrete values on a I?-dimensional torus that we 
refer to as the replica space. The voltages are discretized 
by setting 9 = A9k, where Ad = 9m /M is the gap be- 
tween successive voltages, 9m is a voltage cutoff, fc is a 
D-dimensional integer, and M a positive integer. The 
components of k are restricted to —M < k'^"-^ < M and 
periodic boundary conditions are realized by equating 
^(q) _ fc'-"^mod(2M). The continuum may be restored 
by taking 9m oo and A9 0. By setting M = m? , 
9m = 9q m, and, respectively, A9 — 9o/m, the two limits 
can be taken simultaneously via m — > oo. Note that the 
limit D —> has to be taken before to — » oo in order to 
ensure lim£)_o(2A/)^-^ = 1. Since the voltages and A are 
conjugated variables, A is affected by the discretization 
as well: 



A = AA/ , AAA6I n/M , 



(2.10) 



mod(2MAA) 



and 



{2My 



) mod(2MAe) 



(2.11a) 



(2.11b) 



do hold. Equation ( 2.11 ) provides us with a Fourier 
transform between the 9- and A-tori. It is important to 
note that the replica space Fourier transform of '>Jj^{x), 



^g^{x) = (2M)-^Ee^pO^-^0'^A(2^) 



-D 



(2.12) 



satisfies the condition J^ff^ffi^) — ^ hence is nothing 
more than a Potts spin with q = (2A/)^ states. 

In passing we emphasize the benefit of the replication 
procedure. It provides us with an extra parameter D 
that we may tune to zero. In this replica limit the nor- 
malization denominator Z~^ goes to one and hence does 
not depend on the distribution of the bond conductances 
anymore. Then the only remaining dependence on this 
distribution rests in the power P appearing in the ex- 
ponential in Eq. (2^). In the replica limit, therefore, we 
just have to average this exponential instead of the entire 
right hand side of Eq. (2/7). This average then provides 
us with an effective power or Hamiltonian which serves 
as vantage point for all further calcula tions . The effective 
Hamiltonian will be discussed in Sec. II B. 



Now we come back to the role of Eq. (2^) as a gener- 
ating function. Since the integrations are Gaussian they 
are readily carried out with the result 



G 



(^x, x' , )Cj — P {x, x') / exp 



A2 



R{x,x') 



(2.13) 



It is evident from Eq. (^.13 ) that G{x, a;'. A) represents a 
generating function for the average resistance. To obtain 
a generating function for the average conductance one 
simply needs to carry out a Fourier transformation in 
replica space, 

g(^x,x\9^ = P{x,x')-^J2''''P{'^-^) 



exp 



(2M) 
A^ 



R{x.x') 



(2.14) 



where I is a D-dimensional integer taking the same values 
as fc. This choice guarantees that the completeness and 



After paying due attention to the exclusion of A = C) 
we may approximate the summation in Eq. ( p. 14 ) by an 
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integration, 

g(x, x', 0) (x, x') 



d X exp 



{2MAX)^ 
A2 



R{x,x') + iX 



{2MY 



(2.15) 



The A-intcgration is straightforward since it is Gaussian. 
In the Umit 13 we obtain 



G 



{x,x',6^ =P{x,x') 



exp 



P{x,x') 



-R{x,x') ^ 



-Afl^i {x,x') + 



(2.16) 



We learn from Eq. (2. If) that G{x,x',9) is indeed the 



generating function we are looking for and that Mp^i can 



be extracted simply by taking the appropriate derivative, 

d 



M|?i {x,x') ^P{x,x') 



-G 



(^x,x',e^ 



(2.17) 



We conclude Sec. |II A| by addressing the physical con- 
tent of ip^ix) and its replica space Fourier transform 
^g{x). A reasoning similar to that for the two-point cor- 
relation function G{x, x' , X) leads to 



rep 



Pr. 



exp 



(2.18) 



where the prime now indicates averaging subject to the 
condition that x is located on an infinite cluster. Poo 
stands for the percolation probability that a point be- 
longs to an infinite cluster and Rooix) de notes the resis- 
tance between x and infinity. From Eq. (2.18) we learn 
an important feature of ip^lx), namely that its average is 
proportional to the percolation order parameter Poo ■ For 
reasons that are clear by now it is preferable to consider 
$g-(a;). Upon Fourier transformation in replica space we 
find in the limit D ^ 



exp 



- 1 



(2.19) 



Anticipating results we will derive in Sec. IV we rewrite 
Eq. (pTTgl) as 



P^ 



X exp 



dtp{t) 



32 



1 , 



(2.20) 



where w is a constant proportional to CTq ^ and where 

p{t)^(6(t-we""^^~>'''R^))\ (2.21) 

is the probability distribution of the conductance to in- 
finity. Thus, the physical meaning of the averaged ^g{x) 
may be stated as follows: {<^g{x))^ci> corresponds to the 
percolation order parameter times a scaling function that 
incorporates the distribution of the conductance to infin- 
ity. 

B. Field theoretic Hamiltonian 

This section presents our derivation of a field theoretic 
Hamiltonian for the RRN^*^. It is guided be the work of 
LT. 



We start by revisiting Eq. (2/7) from which we read off 
the effective Hamiltonian announced in Sec. HA: 



In ( exp 



In I y ]^ dcTfc fl (cTfe) exp 



(2.22) 

For the subsequent steps it is convenient to recast 
Eq. ( I2.22D as 



(2.23) 



Here, we have introduced 



K 



lull 



V j da f (cr) exp 
'0 



-ere 



(2.24) 

with V = p/{l — p). Moreover, we dropped a constant 
Nb ln(l — p) with Nb being the number of bonds in the 
undiluted lattice. In order to refine -ffrcp towards a field 
theoretic Hamiltonian we now expand K{0) in terms of 

1 



iX- 



{2M) 



(2.25) 



where K{X) is the replica space Fourier transform of 
K{9). To evaluate K{X) we approximate the summation 
over 9 by an integration. This gives, up to a multiplica- 
tive factor that goes to one for D ^ 0, 



K X\ = 



61 exp 



— 00 







iX-9 


ln| 



/■Co 

/ dcr f (cr) exp 
Jo 



(2.26) 
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Upon expanding the logarithm and carrying out the 9- 
integration we obtain, once more by dropping a multi- 
plicative factor that goes to one in the replica limit, 



K A 



E 

1=1 

with Fi{X) being given by 



(2.27) 



F, 





f 


) ^ r 


dai-- ■ / 




Jo 




A2 


X exp 






2 (Tl 4 



dai f (ai) ■ ■ ■ f (ai) 



(2.28) 



Integration yields, as demonstrated in Appendix 



Fi (a) = 1 + Ai S? + Bi A2'(i-'') + 0[{}?Y 



(2.29) 



where A^"^^ is understood as (A^)'^^ Ai and Bi 
are constants. For example, Ai is given by Ai = [1 — 
a)/{2aao) and Bi reads Bi = -r(a)/(2cro)^"" with T 
denoting the F-function. The general form of the Ai is 
Ai ~ + l{a — 1)]"'^- For reasons that will be given 

in Sec. [I C we will neglect all ter ms a ssociated with i ?;>i 
from now on. By inserting Eq. ( 2.29| ) into Eq. ( 2.27 ) we 
find 



Kit] =K + w}? + t;A2'(i-") + 



(2.30) 



with if, w, and v being expansion coefficients, v is pro- 
portional to ctq"^ ^'^^ positive for a > 0. w is propor- 
tional to (Tq^ ■ Its sign depends on the values of a and 
V. From now on we omit factors (2M)~^ that go to 
one in the replica limit. Moreover, we define the discrete 
gradient V^- via 



(2.31) 



Collecting we find that 



S'(2M)° symmetry is lost in favor of an 0{D) rotational 
symmetry in replica space. 

We proceed with the usual coarse graining step and 
replace the Potts spins ^g{x) by order parameter fields 
(y9(x, 0) that inherit the constraint = We 

model the corresponding field theoretic Hamiltonian Ti. 
in the spirit of Landau as a mesoscopic free energy. The 
constituting elements are local monomials of the order 
parameter field and its gradients in real and replica space. 
Purely local terms in replica space have to respect the full 
^{2M)° Potts symmetry. After these remarks we write 
down the Landau-Ginzburg- Wilson type Hamiltonian 



e 

- |(p(x,0y|, (2.34a) 



with the kernel being given by 



K (A, A,-) ^T-W^ -w^}+v[ - V| 



(2.34b) 



In Eq. ( ^.34| ) we have neglected all higher order terms 
that are irrelevant in the renormalization group sense, 
w, and V are now coarse grained analogues of the orig- 
inal coefficients appearing in Eq. (2.32). The parame- 
ter T — Tc {pc — p) specifies the deviation of the oc- 
cupation probability p from the critical probability Pc- 
In mean field theory the percolation transition happens 
ai T — Tc = 0. We point out that TL reduces to the 
usual field theoretic Hamiltonian for the {2M)^ states 
Potts model upon setting w = v = 0. Thus, Ti. satisfies 
an important consistency requirement since one retrieves 
purely geometrical percolation in the limit oo. 



C. Scaling analysis in the voltage variable 



In the limit of perfect transport, ctq 
cients w and v vanish and iJrcp reduces to 



Now we address the relevance of the coupling constants 
K — i(;V|+ V (^^l) ^g{j) ■ associated with the bond conductances. We carry out a 

scaling analysis by rescaling the voltage variable: 9 b9. 
(2 32) Here b denotes a scaling factor and should not be con- 
fused with the index labeling the bonds. By substituting 
oo, the coeffi- (^(x, 9) — ip'ix, b9) into the Hamiltonian we get 



(2.33) 



<»,i> 



This Hamiltonian represents nothing more than the 
(2M)^ states Potts model that is invariant against all 
(2M) ! permutations of the spin states. If Cf^^ ^ 0, this 



j d^'x J2 { 1^' (x, b9) K (A, A,-) ^' (x, b9) 

9 

f/(.,«) 



(2.35) 
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Renaming the scaled voltage variables 0' = hO leads to 



I d^^xY, (x, e') K (A, fe^A.) ^' (x, e') 



I. J 

6 



(2.36) 



Obviously, a scaling of the voltage variable results in a 
scaling of the voltage cutoff, 9q b0o. However, by 
taking the limit Z? — > and then m — s- oo, the dependence 
of the theory on the cutoff drops out. In other words: Bq 
is a redundant scaling variable. Thus, one can identify 6' 
and 6 and conclude that 



H 



■,T,W,V 



ip{^^,dy,Tyw,b'^^^-''^v . (2.37) 



Next we consider the consequences of Eq. (PItI) for 
the correlation functions of the field <^(x, 0) given by 

Gat ({x-^} ;t,w,v^ 

= j Dip ip{yii,e-i^---ip(^N,dN 

X exp (y-'H ip ^x, (tj ; r. 



(2.38) 



where Vip indicates an integration over the set of vari- 



ables {((^(x, 0)} for all x and 6. Equation (2.37) implies 
that 

Gat ({x, ■,t,w,v^ 

= Gn ({x, 6^} ; r, h^w, b^^^-'^'^v^ . (2.39) 
From Eq. (2.3E) in conjunction with Eq. ( 2.16| ) we deduce 



^Ml^i ((x,x') ;r,ii;,w) 



((x,x');t, 62u;,62(1"'^)w 



(2.40) 



The freedom of choice with respect to b has not been 
exploited yet. To address the issue of relevance we choose 

62 = ^-l/(l-a)^ rpj^ig g-ygg 

M|C, ((x,x');T,w,t;) 

= t;-i/(i-)AfSC^ ((x,x');t,W«'/(i-''),i) . 

(2.41) 



By virtue of w ~ ctq '•^ we may recast Eq. (2.41) as 



Mft, ((x,x');r,w,t;) 



= ao-i/i((x,x');r,«;/t,i/(i-)) , (2.42) 



with /i being a scaling function. We learn that w ap- 
pears only in the combination w/v^^''^~°'\ A trivial con- 



sequence of the fact that the Hamiltonian (2.34) must be 
dimensionless is that wX^ ~ wA^^^""-' ~ where /i is 
an inverse length scale. Thus, w /v'^/'^^~°'^ ~ ^-2a/(i-a)^ 
This leads to the conclusion that w is marginal for a of 
order e whereas it is clearl y irr elevant for a of order one. 



As hypothesized in Sec. II B the scaling analysis in the 



voltage variable j ustifie s that we have neglected in the 
remainder of Eq. ( 2.29|) all terms associated with Biyi. 
Suppose that we had retained these terms. Each of them 
had contributed a term — «;(— V?.)'*^^~"^ to the kernel in 
Eq. ( 2.34| ). From the preceding paragraph it is evident, 
however, that vi appears in the average conductance only 
as vi/v'- ^ We conclude that keeping the i?i>i had 

produced only irrelevant terms and that neglecting them 
in studying the leading behavior at the critical point is 
indeed justified. 

For our RG improved perturbation calculation pre- 
sented in Sec. IV it will be helpful to dispose of a coupling 
that is invariant under —>■ bO. To identify a candidate 
we revisit Eq. (2.40) and choose 6^ = w~^. This leads to 



Mf^i ((x, x') ; T, w, v) = (Tq- V2 ((x, x') ; r, /i) , 



(2.43) 



with /2 being another scaling function, h = v/w^^"" ~ 
^2a ^m-Qg out to be the sought invariant coupling con- 
stant. We will see that it emerges quite naturally in 
perturbation theory. Hence, we refer to h as effective 
coupling. 



III. REVIEW OF THE RRN 

This section he re p resents a brief review of the model 
with a = [|l|, |lg |0|. We provide the reader with 
background on the RRN to make the subsequent analysis 
of the RRN^^ more digestible. 

In this as well as in the following section we utilize Ti 
and calculate the generating function G(x, x',6') by em- 
ploying field theory augmented by renormalization. For 
background on these methods we refer to Ref. As 
soon as we have G(x, x', 9) at hand it is a straightforward 
matter to extract M^^. 

For a — the coupling constants v and h are redundant 
and can be set to zero. Straightforward dimensional anal- 
ysis shows that dc = 6 is the upper critical dimension and 
that r2 and Fa are the only superficially divergent vertex 
functions. The diagrammatic elements as constituents of 
our perturbation calculation are the three-leg vertex g 
and the principal propagator 



G^°'^(p,A) = G(k,A){l-(5^g}, 



(3.1) 



where G(p, A) — (r -I- + wX^) ^. Due to the fac- 
tor {1 — 6^^} that enforces the constraint A ^ the 
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principal propagator decomposes in a conducting part 
G"^°"'^(p,A) = G{p,X) carrying replica currents and an 
insulating part G™^{p) = G(p, X)d^ g not carrying replica 
currents. Each principal diagram decomposes into a sum 
of conducting diagrams consisting of conducting and in- 
sulating propagators. 

Our real- world interpretation H, |l9|, |2[ ||, H, |2|, 
p6| , |27[ |, in which the conducting diagrams are viewed as 
being resistor networks themselves, provides for a pow- 
erful and elegant framework to calculate these diagrams. 
At first we rewrite the propagators in Schwinger param- 
eterization, 



G'(p,A) 



ds exp 



(3.2) 



Next we interpret the Schwinger parameters s of the con- 
ducting propagators as their resistance. Then we can 
express the A-dependent part of any conducting diagram 
with P'^""^ conducting propagators in terms of its electric 
power P: 



exp 



^ SiXf j = exp 



wP 



(a, {k} 



(3.3) 



The summation on the left hand side runs over all con- 
ducting propagators. A^ = Ai(A, {k}), where A is an ex- 
ternal current and {/?} is a complete set of independent 
loop currents, denotes the current flowing through con- 
ducting propagator i. In this representation it is easy to 
see that the sum over the loop currents is determined by 
the total resistance R{{si}) of the respective diagram, 



E 



exp 



wP 



(a, {4) 



exp 



-Ri{s.,})wX^ 



(3.4) 



Carrying out the usual momentum integrations, which is 
straightforward after completion of squares, and Taylor 
expansion gives 

/ (p2, A^) = Ip (p2) - (p2) wX^ + --- 



Y[ds, [l-i?({s,})u;A'-|- 

i 



(3.5) 



for the overall form of any conducting diagram. 
D(p^,{si}) stands for the usual (Schwinger parameter- 
ized) integrand of the corresponding diagram in the stan- 
dard if^ theory. 

The ultraviolet (UV) divergences encountered in com- 
puting the diagrams can be handled by resorting to di- 
mensional regularization. In dimensional regularization 
the UV divergences appear as poles in the deviation 
e = 6 — d from dc- These poles may be eliminated from 



the superficially divergent vertex functions by employing 
the renormalization scheme 



r 



9^9 



= ^-3/2^1/2^71/2^1/2^^/2 



(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 



where the ° indicates unrenormalized quantities. The am- 
plitude Ge = (47r)~'^/^r(l -|- e/2) is introduced for con- 
venience. Z, Zr, and Zu are the usual Potts model Z 
factors known to three-loop order . 

In the minimal renormalization procedure, i.e., dimen- 
sional regularization in conjunction with minimal sub- 
traction, the Z-factors are of the form 



Z..Xu) = 1 + 



E 

m— 1 



x(")(i 



(3.7) 



The X ('")(«) are expansions in the coupling constant u 
beginning with the power u™. It is a fundamental fact of 
renormalization theory, cf. Ref . , that this procedure 
is suitable to eliminate the UV-divergencies from any ver- 
tex function order by order in perturbation theory. 

The unrenormalized theory has to be independent of 
the arbitrary length scale introduced by renormal- 
ization. Hence, the unrenormalized correlation functions 
satisfy the identity 



{x, a| ;5,f. 



w 



0. 



Equation (3^) translates via the Wilson functions 



7... (u) 

K (u) 



^ 1 7 

In A. 

du 
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9/i 
dlnw 



u (37 - 7„ - e) 



7-7r 



1 -In 



(3.8) 

(3.9a) 
(3.9b) 
(3.9c) 
(3.9d) 



(the |o indicates that bare quantities are kept fixed while 
taking the derivatives) into the RGE 



d_ 

X Gat 



■(3 



d 



du 



d 

OT 

u, r, w, fi 



dw 

: 0. 



N 
+ ^7 



(3.10) 



Fro m the structure of the renormalization factors 
( |3.13D and the Wilson functions one can deduce 

the important fact the RGE is determined entirely by 
the X^^Xu)- This may be seen as follows. From the 



definitions (3.9a) and (3.9b) we learn that the Wilson 7 
functions can be expressed as 



7... (m) 



P in) l^Z.. 



(3.11) 
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A glance at Eq. (3.7) tells us then that the logarithmic 
derivative in Eq. (3.11) has a pure Laurent expansion 
with respect to e starting at first order in 1/e. Moreover, 
we know from Eq. ( 3.91 ) that begins with the zero- 
loop term —eu. Because the 7 functions are finite for 
£ ^ their e poles have to cancel order by order in the 
loop expansion. As a consequence, the 7 functions are 
given by 



7... H 



(3.12) 



For this reason we will focus in the remainder of this 
paper on the X^^u). We neglect higher order terms in 
the expansion (3.7) and write 



00 y(L) 

Z,„{u) = 1 + ^ u^ + O (s-^^ 



L = l 



Le 



(3.13) 



where the y(^) are numerical coefficients independent of 
e. 

The RGE can be solved in terms of a single flow pa- 
rameter I by using the characteristics 



dl 








(3.14a) 


du 
~dl 




u(l) 


= u, 


(3.14b) 






f(l) 




(3.14c) 


/— Inw 
dl 


= C("(0) , 


w{l) 


= W , 


(3.14d) 


'Jl'-' 


= lim , 


Z{1) 


= 1. 


(3.14e) 



These characteristics describe how the parameters trans- 
form if we change the momentum scale fi according to 
fi — > p,{l) = Ifj,. Being interested in the infrared (IR) 
behavior of the theory, we study the limit 1^0. Ac- 
cording to Eq. ( 3.14h| ) we expect that in this IR limit 
the coupling constant u{l) flows to a stable fixed point 
u* satisfying (3{u*) — 0. The IR stable fixed point so- 
lution to the RGE is readily found. In conjunction with 
dimensional analysis it gives 

G^({x,A};u,r,z«,^) 

X Gat (|/x, a} ; u^'J-^/^t, r't'/'^w, /i) 

(3.15) 

with the critical exponents for percolation rj = j(u*) and 
v = [2 — k{u*)]~^ known to third order in e |28| . — 
i>[2 — Ciu*)] is the percolation resistance exponent known 
to second ord er in e p^ . 

Equation ( [3.15| ) implies that the two-point function 
G(x, x',0) scales at criticality r = as 



G 



(x, x', = i^'^/^n |x - x'l , r^/'^w-'e') 



(3.16) 



with (3 = {d — 2 + r])i'/2 denoting the percolation or- 
der parameter exponent and where f2 is a scaling func- 
tion. Upon choosing / — |x — x'|~^ and expanding the 



right hand side of Eq. ( 3.16 ) we obtain by dropping non- 
universal constants 



G(x,x',0^) =|x-xf''/'^ 
x< -^„-l||x-xr*/''- 



(3.17) 



With help of Eq. (|t^ it is now straightforward to de- 
duce the scaling behavior of the average conductance. 



w 



(3.18) 



IV. 



RENORMALIZATION GROUP ANALYSIS 
OF THE RRN^° 



Now we turn to the RRN^^ and assume that < a < 
1. First we address the renormalization of the model. 
By carefully analyzing the RG flow we reveal the critical 
behavior of and S^^. As far as notation is con- 

cerned, we adopt the same convention as in Sees. I and 
||. Quantities that might be confused with their analogs 
for the RRN are marked by the superscript SC. 



A. Renormalization 

Obviously the Potts model Z factors are independent 
of a and hence do not require further consideration. The 
Z factor pertaining to w, however, is expected to be dif- 
ferent from its analog for the RRN, i.e.. 



w 



w 



z-'^z^^ 



w . 



(4.1) 



Since u 7^ we need an additional renormalization 
v^v = Z-^Z^'^v. 



(4.2) 



From Eqs. (4.1) and (4.2) we deduce immediately, that 
the effective invariant coupling h — v/w^~'^ announced 
in Sec. II C has to be rcnormalized by 



h = z-'^zf^{zl^X-'h^^ 



(4.3) 



The factor is included to render the rcnormalized 
effective coupling dimensionless. 

After these remarks we now bring our attention to the 
perturbation calculation and its Feynman diagrams. For 
the RRN^^ the principal propagator has a form similar 
to its analog for the RRN, 



^SC,bold(p^X) = GSC(k,A){l-%g}, 



(4.4) 



However, its flesh is now given by G^*-"(p, A) = (t -I- -I- 
wA2-|-t;A2(i-''))-i. Evidently, GS^'^°'^(p, A) decomposes 
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into a conducting and an insulating part. This leads to 
conducting diagrams identical to those for the RRN up to 
apparent distinctions in the definition of the propagators. 
Due to these distinctions an expansion for small external 
momenta and currents leads to 



Iv vX 



■2(l-a) 



(4.5) 



instead of Eq. ( |3.5[ ) for the overall form of the conducting 
diagrams. 

To extract information on Z^'~^ consider the A depen- 
dent part of a conducting diagram with P'=°"'^ conducting 
propagators which reads in Schwinger parameterization 



{ft} 



exp 



wXt + vX, 



2(l-a) 



(4.6) 



We keep in mind that Xi — Xi{X, {/?}). The important ob- 
servation is now that any conducting propagator affected 
by the summation over the loop currents gives a contri- 
bution to Eq. (4.5) polynomial in A, i.e., it contributes 
to Iw rather than to ly- The only contributions to ly 
can come from conducting propagators not affected by 
the summation over {k}. In the terminology of our real- 
world interpretation this means that ly is determined 
exclusively by the red bonds of that diagram, i.e., by its 
singly connected conducting propagators. To be specific, 
Iv is given by 



pQC 



(4.7) 



where the sum runs over all P'^°'^ conducting propagators 
of the diagram not belonging to any closed conducting 
loop. 

Now we take a short detour and recall some central fea- 
tures of our field theory on the nonlinear RRN 
in which the occupied bonds obey a generalized Ohm's 
law V ^ I^. The field theoretic Hamiltonian for the 
nonlinear RRN corresponds to that for the standard 
RRN with u;V| replaced by Wr Ef=i(-3/56'("))''+i. Ac- 
cordingly, one encounters a generalized renormalization 
factor Ziu^ that then leads to a generalized resistance 
exponent (j)r — v\l — Crju*)] . Th e Wilson function 
is defined analogous to Eq. ( [3.9d ) with 7^, replaced by 
Iwr — ln.^M)rlo- The generalized resistance exponent 
has the physical meaning of governing the average non- 
linear resistance at criticality. 



Mr{x,x') = {x{x,x')Rr{x,x'))c / {X{x,x'))c 



(4.8) 



The nonlinear RRN is particularly interesting, because 
it provides for an elegant way to determine fractal dimen- 
sion of percolation clusters by considering specific values 
of r. For example, it is a well known fact |3l|] that 



hm Mr 



M, 



red ; 



(4.9) 



where Mrod is the mass (aver age number) of the red 
bonds. From Eqs. ( |4.S| ) and (O) one obtains imme- 
diately the scaling relation dred — limr_+oo '^rlv for the 
fractal dimension d^^A of the red bonds. It was shown 
rigorously by Coniglio js^ that dred = ^lv which means 
(Til,. = 1. From the definition of Sr it then 



that lim,.^oo 0r 
follows that lim,. 
the identity 



, Cr(u) = k(u) which leads in turn to 



lim Zu,^ — Zt- . 



(4.10) 



In Refs. |2^, ^ we showed based on our real- world in- 
terpretation that the contribution Iw^ of a conducting 
diagram to Z^^ takes for r — > 00 the simple form 



/•oo 

lim%, (p2)= / TTds, V s,i?(p2,{s,}) 



Comparison of Eqs. (4.7) and (4.11) yields 
lim Z,„„ = . 



(4.11) 



(4.12) 



From this and Eq. ( 4.10| ) we finally conclude the identity 

(4.13) 



Zf^ = Zr . 



For analyzing Z'^ we revisit Eq. (4.6) and rescale the 
replica currents, A^ — ^ ui^^A^. This recasts Eq. (4.6) into 



E 

{«} 



exp 



E K 



Xi -f h\ 



■2(l-a) 



(4.14) 



As mentioned above, conducting propagators belonging 
to closed conducting loops (blobs) lead in an expansion 
for small A to terms polynomi al in A that contribute to 
Iw ■ Since h appears in Eq. (4.14) /^y will in general 
depend on h (cf . Appendix ^ . We conclude that Z^ is 
not only a function of u but also of ft,. 



Zl'' ^Zl^(u,h). 



(4.15) 



For arbitrary a G (0, 1) it is difficult to gain further 
insight into Z'^ . Anyway, the most exciting values of a 
are those for which a is of the order of our small expansion 
parameter e. This is the key region in which the crossover 
between the RRN and the RRN^^ occurs (the naive limit 
a ^ presupposes e <C a and hence is inadequate to 
resolve the crossover). In case both a and e are small 
we can follow the work of Honkonen and Nalimov |33j to 
analyze the structur e of Z ^p . This structure differs from 
that of Zw [cf. Eq. ( ^.13 )] because two major differences 
emerge for a > 0. First, the conducting propagators 



°° \2(l-a)fe 



t=0 [t + p2 + A2] ■ 



(4.16) 
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give rise to an entire series of individual terms. Second, 
since a is now of order e, the poles in these individual 
contributions are now of the type l/(i£ + 2ka). Hence 
Z^p is of the form 



Obviously value of ^ at the fixed point u* is given by 



{L,k) 



l,k 



L=l k=0 



Le + 2kc 



■u'h 



O e- 



(4.17) 



Here, the yi^''^^ are the numerical coefficients of the 
poles. Appendix ^ illustrates the preceding arguments 
at the instance of an one-loop calculation. 

In minimal subtraction, the Yw^'''^ are independent of 
both s and a. This fact provides us with a simple method 
to calculate the numerical coefficients because it is suf- 
ficient to consider the limit a ^ 0. In this limit the 
Hamiltonian for the RRN reduces to the Hamiltonian 
for the RRN with w replaced by w + v. As a consequence 
we obtain by exploiting Eq. ( 4.13 ) the relation 



Zl^iu, h) = Z.^{u) + h [Z.^{u) + Zr{u)] (4.18) 

valid i n th e limit a 0. I nserti ng the explicit 
forms ( ^.131) and (|4.17D into Eq. ( [4.18]) we obtain the 
following relations between the numerical coefRcients: 

^(L,0) _ ^(L) ^(L,l) _ ^(L) _ ^(L) , ^(L,k>l) _ „ 

Recalling that the Fi^''^^ are independent of a we con- 
clude that 



zSC(u,/,) = 1 + ^^ 



Y 



y^(L) ^(L) 



L=l 



Le 



Le + 2a 



(4.19) 



B. Scaling 



We proceed in the same fashion as in Sec. HI and set 



up a Gell-Mann-Low RG equation for the RRN^^ . By 
carefully analyzing the RG flow we derive the scaling be- 
havior of the average conductance. 

The RGE for the RRN^"^ is somewhat richer than that 
for the RRN: 



— +f3— 
^ djjL du 



d 



where we have introduced the Wilson functions 



91ni 



dfj, 
dlnv 



N 

^\ 
(4.20) 



(4.21a) 
(4.21b) 
(4.21c) 



\lv. 



(4.22) 



In order to express at u* in terms of the known RRN 
exponents we introduce the function 



which leads to 

cSC(ii*,;j) = 2-0/i.-/(/i) 



(4.23) 



(4.24) 



Equations ( [4.2l| ), ( [4.22| ), and ( [4.24| ) provide us with flow 
equations for the couplings w and v: 



4>/,^~f{h{i)), w{^) 



/|ln.(0 



2-1/v, v{l) = v. 
The flow equation ( 4.25b| ) for v is readily solved, 



v{l) = v/^-i/^ 



w, 

(4.25a) 
(4.25b) 

(4.26) 



To solve the flow equation ( 4.25a ) for w we first have to 
analyze the flow of h. 



From the renormalization of h (4.3) follows immedi- 
ately that the logarithmic derivative of the renormalized 
h with respect to fj, is given at u* by 



i—u* 1 
V 



+(l-a)/(/i). 



(4.27) 



Consequently, h obeys the flow equation 

9.... .... n 



h{l)\ - (l-a)0-l 



1 



^)/(M0) 

(4.28) 



with the initial condition ft,(l) = h. A glance reveals that 
the flow of h has two IR fixed points, viz. h\—Q and h\ 
determined by 



1 



1 - a 



(4.29) 



Now that we know the fixed points of h we have to 
analyze their stability. Since there are two fixed points 
either one will be stable and the other one will be unsta- 
ble depending on the value of a. Consider h\. h\ is stable 
if the {• • • } bracket on the right hand side of Eq. ( 4.28| ) 
has a positive sign and it is unstable if this bracket has a 
negative sign. For a of order one the sign is certainly neg- 
ative. This leads to the conclusion that is stable for a 
of order one. In the crossover region, i.e., for a of order e, 
the question of stability is more intricate. Here, we need 
more information on the functional dependence of / on 



11 



h. This informatio n can be derived from the st ructu re of 
Z^^ given in Eg. (H.19| ). Upon inserting Eq. ( [4.19[ ) into 
the definition (4.21a) of 7^^ we obtain 



LXJ 



L=l 



y(L) _ y(L) 



(4.30) 



Substituting ^^'^{u*,h) into the definition (|4.23| ) of / 
leads to 



/(/i) = 1] 



(4.31) 



The important conclusion from Eq. (4.31) is that / is 
linear in h for a of order e. Hence, it r esult s in a con- 
tribution to the right hand side of Eq. ( 4.28 ) quadratic 
in h that can be neglected in analyzing the stability of 
h^. Finally, we perceive that /ij is stable ii (p > (1 — a)~^ 
whereas /12 is stable if0<(l — a)^^. 

At this point we possess enough information to solve 
the flow equation (4.25a) for w. For > (1 — a)^^ we 



know that / tends to zero because h flows to /i*. For 
< (1 — o)~^, on the other hand, we have to insert 
(4.29 ). Sum marizing both cases we write the solution to 
Eq. ( |4.25aD as 



(4.32) 



where we have defined the a dependent resistance expo- 
nent 



if 



if 

1 — a 



> 
< 



1-a 
1 



(4.33) 



Collecting t he ab ove results we find that the solution 
to the RGE ( 4.20 ) augmented by dimensional analysis 
reads 

GN[{^,x};u,T,w,v,fi) =1'""'' 

X Gat ({/x, a} ; w*, r^/^^r, /"'^"''('^^/'^w, r^'^v, ^i) . 

(4.34) 

This scaling form implies for the Fourier transformed 
two-point function at criticality that 

xs(?|x-x'|,r'^''^(")/''ii;-i0^,;it,2) , (4.35) 

where S is a scaling function. The choice I = |x — x'|^^ 
and subsequent Taylor expansion yields 



G 



(x,x',0^) 



= |x — x 

1*2 



,,2/3/1/ 



X < - w — X - X 



(4.36) 



With help of Eq. ( ^.17 ) we readily obtain 



(4.37) 



It remains to deduce the scaling behavior the average 
conductivity and its conductivity exponent t^'~'{a). Com- 
monly, the conductivity of percolating systems is defined 
with respect to a bus bar geometry where the network is 
placed between two parallel superconducting plates (the 
electrodes) of area L'^~^ a distance L apart. From the 
above we expect that the average conductance of 
this system scales as 



(4.38) 



where Ha is an a-dependent scaling function with the 
properties 



const for x 1 
f-Qj, x > 1 



(4.39) 



Now consider length scales large compared to the cor- 
relation length ^. In this regime the RRN above the 
percolation threshold, r < 0, can be viewed as a homo- 
geneous system of conductivity S^^(r). Hence, we may 
write for L 3> ^ that 



I]SC(r)^L2-VC(L,r) 



By virtue of Eq. ( 4.39 ) we finally get 



sSC(r) ^ |r 



(d-2)iy+0^*^(a) 



(4.40) 



(4.41) 



which means that the conductivity exponent is given by 
t^^{a) = id-2)iy + (f>^^ia). (4.42) 

V. CONCLUDING REMARKS 

We showed without relying on the assumptions of the 
'nodes, links, and blobs' picture that t he c onductivity 
exponent for the RRN^^ is given by Eq. (L5). For suffi- 
ciently large values of a this means that 



(5.1) 



recognized by Halperin et al. as a lower bound to 
tS^(a) that can be attributed to a dominance of the red 
bonds. At a critical value = 1 — 1/0 the conductivity 
exponent crosses over to the value 



(5.2) 



for the standard lattice model and is essentially deter- 
mined by the blobs. 

The 'nodes, links, and blobs' picture provides an in- 
tuitive explan ation for this behavior. Since the distribu- 
tion function (2.1b) is singular, there exists at criticality 
a large number of bonds with arbitrarily large resistance. 
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Whenever one of these bonds is red it may dominate the 
total resistance of its Hnk. If the large resistance, on 
the other hand, occurs in a blob it has, in general, lit- 
tle impact because the current can flow through parallel 
paths. For a > Qc the resistance of the weakest red bond 
is typically larger than the sum of resistances of the other 
building blocks (the blobs and the other red bonds) of a 
link so that the weakest red bond dominates the total 
resistance of its entire link. For a < ac the importance 
of the weak red bonds is diminished and the network be- 
haves effectively as a standard RRN. 

We point out that our analysis was not restricted to 
any particular order in e-expansion. Unlike LT who set 
up a 'momentum shell' RG, we used the powerful meth- 
ods of renormalized field theory. These methods allowed 
us to explore general properties of the renormalization 
factors, subsequently the RGE, and finally the conduc- 
tivity exponent to all orders in perturbation theory. 



Acknowledgments 

We acknowledge the support by the Sonderforschungs- 
bereich 237 "Unordnung und grol3e Fluktuationen" of the 
Deutsche Forschungsgemeinschaft. 



APPENDIX A: CONTENTS OF THE COUPLING 
CONSTANTS w AND v 

In this Appendix we give details on the integrations in 
Eq. ( 2.2§| ). At the instance of Fi and F2 we illustrate the 
contents of the coupling constants w and v. 

We start by recalling the definition of Fi. After 
changing the integration variable from cr to t by setting 
t = A^/(2(t) Fi takes on the form 



1-a 




Here 



r(a,y) 



(Al) 



(A2) 



is the incomplete Gamma function that has about y = 
the Taylor expansion |34 



r(a,y)^r(a)-^ ^ ,y I . 

nl (a + n) 



(A3) 



Now we turn to F^. After isolating the contribution for 
A = it is save to change variables by setting u\ = u^tx 
and (72 = (70^(1 — a;). We get 

F2(y^ = l + (l-a)2^ dtj^ dxt^-^^x-^il-x)-" 

^{^^p(-2S)-4- ^^'^ 

The integration over x gives a Beta function |Q B{1 — 
a, I — a). The integration over t can be handled analogous 
to the integration in the preceding paragraph. As a result 
of both integrations we obtain 



F, 



(X^ = l + {l-afB{l-a,l-a)l - 



22(l-a) 

2(1 -a) 



+ 



2an 



2(l-a) 



A2 



With help of the expansion ( [A^ ) we finally get 



F2IX) = l + {l~ayB{l-a,l-a)i - 



2(1 - a) 



X r(2(a-l)) 



4cto 



2(l-a) 



1 A^ 

2(a - 1) Aao 

(A7) 



The higher orders i^;>2 can be analyzed by similar 
means. It is not difficult to convince oneself that the 




O A 



general form of the Fi is given by Eq. (|2.29| ) and that 
the coefficients Ai decorating A^ are proportional to 
[1 + l{a — 1)]~^. Since the Fi enter into the kernel K{\) 
via Eq. ( ^.27|) the coefficient of the A^-term in K{X) is 
given by 



1 = 1 



E 

1=1 



Ci 



I l + l{a-l) 



(A8) 



where the C; are po sitiv e constants proportional to . 
We learn from Eq. (A8) that the value of w depends on 
details like the specific values of v and a. It may be finite 
or infinite, positive or negative, v depends on v and a as 
well, see Eqs. (2.27) and (A4). However, for < a < 1 



and p > it is always finite and positive. 



Thus, we obtain 
Fi (a) ^ 



i-r(a) 



A2 \ 1 - a A2 



2(70 



+ 0{{S?f) 



a 2cto 



APPENDIX B: ONE-LOOP CALCULATION 

In this Appendix we calculate the one-loop contribu- 
tions to explicitly. We will see that the one-loop 



(A4) 

result is in conformity with Eq. (|4.17|) 



13 




FIG. 1: The principal one-loop self-energy diagram (bold) 
decomposes into the conducting diagrams A and B assembled 
of conducting (light) and insulating (dashed) propagators. 



In one-loop order there exists only a single principal 
self-energy diagram that decomposes into conducting di- 
agrams as shown in Fig. |^. In the following we set exter- 
nal momenta to zero for simplicity. Hence the one-loop 
contribution p^"'°°P to the vertex function r2 (note that 
the vertex functions Fjv are defined as the negative of the 
corresponding diagrams) is given by 



2B- A. 

Applying the usual Feynman rules yields 



(Bl) 



F 



1-loop 



2 



^G"=°""(q,K)G™""(q,K + A 

1 K 

cond 



q,A G 



(B2) 



where we have dropped the superscript SC because we 
consider exclusively the RRN^*^ throughout the entire 
Appendix. / is an abbreviation for {2'k)~'^ dfiq. For 



the following steps it is convenient to recast Eq. (B2) as 

,2 r ' n2 



F^-'°°p(a 



9_ 

' 2 



E 



q,K 



G™"'*(q,A)G'"'^(q) 



q, K 



(B3) 



The evaluation of the A-independent term is straightfor- 
ward because E« G™"<^(q, k)^ = G'^^q)"^ for D ^ 0. 

The A-dependent terms are expanded in the external cur- 
rents. After some algebra we obtain 

7A2 + t,A2(l-'^) 



-il-loop 



«4 



1 



(r-f q2 



+ 9' 



E 



[w + {1 - a)vK-^'']\x ■ kY 
(r -I- q2 -I- -I- t;K2(l-i)~)'' 



(B4) 



The first two integrations are readily carried out using di- 
mensional regularization. The summation over the loop 
current can be simplified by exploiting the rotational 
symmetry in replica space and by rescaling wk^ 
After these steps we arrive at 



K . 



^l-loop 



> Gf: 



+ g^wX^I + ---, (B5) 



where 



I = — 
D 



(T + q2+K2 + /lK2(l-o))^ 



(B6) 



remains to be evaluated. Upon rewriting / in Schwinger 
parameterization, the momentum integration can be car- 
ried out immediately. The summation over k can be 
approximated by an integration since we have already 
excluded A = properly. Recasting this integration in 
spherical coordinates we get for D 



I 



1 



6 {Any/^ 



ds "^l"^ exp (-st) 



X [l + (l-a)/iK-2'i]2exp -s{k^ + hK^^^-°- 

(B7) 

The integration over radius variable k is simplified by 
introducing an integration variable x = sk^. Moreover, 
we apply the binomial formu la a nd expand the rightmost 
exponential function in Eq. (B7) with the result 



1 r°° 

^1 d..-^/^exp(-.r) 



12 {iwy/ 

xfi:(:)a-») 

X / dx exp(-a;) 

^0 



Ml! 

fc! 



(B8) 



The integral over x is divergent for a > 1/2 signifying 
that the expansion in A2 is justified only for a < 1/2. 
We assume that a is sufficiently small and abbreviate the 
value of this integral by Ck,n- Changing the summation 
index k to k — n we obtain 

poo 

I = "^Ckh^ dss2-'^/2+afcexp(-sT) , (B9) 
k=o ■'° 

where we have set 



1 



Cfc = 



12(47r)''/2 ^ \n 



■E 



{I -ay 



(-1) 



k — n 



(fc-n)! 



Cfe- 



(BIO) 



The integration over s yields 

oc 

/ = Y,Ckh''T(- + ka) r-(^/2+fca) (g^^) 



k=0 



For a of order e we can expand the F function which 
provides us with the final result 



/ = 



E: 

k=0 



Ckh^ 
s + 2ka 



.-(e/2+ka) 



(B12) 



where Ck abbreviates Ck — 2cfeF(l -I- e/2 + ka). We see 
that the UV singularities manifest themselves in a series 
of poles of the from h'^ / {Le + 2ka) with L = 1. In other 
words, our one-loop example is in agreement with the 
general insights presented in Sec. |^ 
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